Further quantum-gate methods using selective displacement of trapped ions 
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We consider quantum gates for trapped ions using state-selective displacement of the ions. We 
generalise earlier work in order to treat arbitrary separations between the traps. This requires the 
impact of anharmonicity arising from the Coulomb interaction to be estimated. We show that its 
effects are always small enough to allow high fidelity. In particular, the method can be applied to 
two ions in the same trap. We also show that gates between non-neighbour ions, and hence a Toffoli 
(three-qubit controUed-NOT) gate, can be achieved. We discuss how the gate can be applied to 
logical qubits encoded in the decoherence-free-subspace {|01) , |10)}, where each pair of ions stores 
a single qubit. We also suggest alternatives to the spin-echo method to suppress unwanted terms in 
the evolution. 

PACS numbers: 03.67.-a, 42.50.-p 



I. INTRODUCTION 

This paper discusses the implementation of quantum 
logic gates in trapped ions, using a mechanism introduced 
by Cirac and ZoUer in Ref. [1], in which the ions are 
displaced by a state-selective force, and acquire state- 
dependent phases from their Coulomb interaction. In a 
previous paper [2] we discussed in detail the realization 
of a two-qubit phase gate 



-P{e) = |00> (00| + |01) (01| + |10> (10| + e'" 111) (11 



(1) 



implemented between two ions confined in separate traps. 
When 6 = TT this is the gate '-'Z which is equivalent to the 
controUed-NOT gate {'"X) up to single-qubit Hadamard 
rotations H, where X and Z are the Pauli ax and az 
matrices. In this paper we extend the scheme of the two- 
qubit phase gate in several ways. 

First, we present some methods relevant to the two- 
qubit phase gate. These are the use of different laser 
frequencies and/or switching force directions to elimi- 
nate undesired rotations, and the application of the gate 
to manipulate logical information encoded in a simple 
decoherence-free subspace. 

Next, we consider the three-qubit phase gate ^'~^Z 
(controUed-controUed-Z gate) whose effect is 



(2) 



where a,(3,j € {0,1}. This gate is equivalent to 
the Toffoli gate '^'^X (controUed-controUed-NOT) in the 
sense that it suffices to add two single-qubit Hadamard 
rotations to obtain the Toffoli gate. The combination 
of '^''^'Z with the Clifford group (generated by Z, H,^ X) 
is a universal set for quantum computation. In this re- 
spect ^^Z is similar to '^P(7r/2). Both and '^P(7r/2) 
are used in networks to achieve fault-tolerant implemen- 
tation of a universal set of gates on qubits encoded in 
quantum error correcting codes [3, 4]. They are there- 
fore desirable features of any quantum computing sys- 



tem. We study the problem of how to use the state- 
selective displacement of trapped ions (the "pushing" 
method) to achieve the transformation in Eq. (2) with 
three ions in adjacent microtraps. Our solution involves a 
standard decomposition of '~^'~^Z into five two-qubit gates 
(two controUed-NOT and three controUed-phase ^P{9)). 
In this decomposition, one of the two-qubit gates is be- 
tween non-neighbouring qubits. The interesting feature 
of the pushing method under discussion is that this gate 
between non-adjacent ions can be achieved without the 
need to rearrange the ions or swap the quantum informa- 
tion between them. 

Finally, we extend the discussion of two-ion gates to 
the case where the ions are much closer together than 
was previously assumed. In this regime the Coulomb in- 
teraction introduces non-negligible anharmonicity in the 
confining potential. We quantify this effect and show how 
the anharmonicity influences the sensitivity of the gate to 
thermal motion of the ions. This more general treatment 
includes the case where two ions are confined in the same 
trap. We find that the effect of anharmonicity remains 
small even in this limit, and therefore the gate remains 
comparatively insensitive to thermal motion. 



II. TWO-QUBIT PHASE GATE 

In this section we briefly discuss the two-qubit phase 
gate considered in Refs. [1, 2, 5], in order to give the 
background to the rest of the paper. 

Let us consider a system of two interacting qubits. 
Suppose the evolution is of the form 



\a(3) 



(3) 



where \aP) are states in the computational basis 
{|00), |01), |10), 111)}. The phases Qap have values aris- 
ing from the internal energy of each qubit, and the in- 
teraction energy. We wish to obtain a two-qubit phase 
gate, which produces an overall phase i? if and only if 



2 



both qubits are in the logical state that is 
|a/3) — > \af3), a/J = 0, 



\ap) 



(4a) 
(4b) 



In order to obtain this we next apply local operations 
(single-qubit rotations) on both qubits 



S = Si ^ S2 , 
where we will assume 

51 = |0)i(0|e'^° + |l)i(l|e^^^ 

52 = 10)2(01 6'^° + 11)2(11 e^^^ 



(5) 

(6a) 
(6b) 



The sequence SG corresponds to the phase gate (4) if the 
following set of algebraic equations is satisfied 



Ao + Bo 
Ao + Bi 
Ai+Bo 



Ooi 



0, 
0, 
0, 



where the overall phase of the gate is given by 
§ = Ai+Bi + Qii. 



(7a) 
(7b) 
(7c) 



(8) 



Eqs. (7) are three independent equations in four indepen- 
dent parameters {Aq, Ai, Bq, Bi). This means we can 
choose one parameter. Let us choose A^ = — 0oo/2. 
Then, using Eqs. (7), we get for the other three parame- 
ters 



Bo 
Ai 



-eoo/2, 
-610 + 600/2, 
-601 + 600/2 , 



and the overall phase reads 

^? = 611 - 610 - 601 + 600 . 



(9a) 
(9b) 
(9c) 



(10) 



When 1? = TT we have the quantum logic gate '-^Z which 
is equivalent to the controllcd-NOT gate '-^X. 

Next we briefly discuss the spin-echo or 7r-pulse method 
which can be used to suppress an unwanted term in the 
evolution and thus improve the fidelity (precision) of the 
gate. This method relies on the assumption that we gen- 
erate the same unwanted term in two successive opera- 
tions G. A pair of tt pulses on qubits 1 and 2 produces 
the effect 



where 



i? = i?i (g) i?2 



Ri = \0)i{l\ + \l}i{0\. 



(11) 



(12) 



Wc replace the original gate sequence SG by S'{RG)^ 
and obtain 



\aP) 



\al3) , a/3 = , 



\aP) , aP=l, 



(13a) 
(13b) 



where i? is defined in Eq. (10) and the single-qubit rota- 
tions S' are 



S' = S[^S'2, 



with 



S' 



+ |l)^(l|e'''/^ 



(14) 



(15) 



where we have dropped a global phase —(601 + 6io)/2. 
Prom Eq. (13), the logical state |11) is rotated during 
the new gate sequence by 2-d, so to achieve the desired 
evolution in Eq. (4) we require 2^ = n. This means that 
in the new sequence S'{RG)'^ each operation G will typi- 
cally (assuming a linear dependence in the gate time) be 
applied for a time half as long as in the original sequence 
SG. 



A. Fidelity 

The calculation of fidelity of a two-qubit phase gate 
was discussed in detail in Refs. [2, 6]. Here we briefly 
present the results, in order to clarify the importance of 
the TT-pulse method, which will be extended to a three-ion 
gate in Sec. IV. 

We consider a general initial state 



1*0) = X] Cc«/3|Q!/3) 



(16) 



a, 13=0 



and calculate the fidelity minimised over 0^/3, i.e. we 
consider a "worst-case" initial state of the qubits. The 
fidelity of a two-qubit phase gate corresponding to the 
pulse sequence SG can be expressed then as 



mm 



E 

a,/3=0 



(17) 



where SOa/s is a random departure of the phase O^^ from 
its expected value (6^/3) and (•) denotes averaging over 
all degrees of freedom of the system. For trapped ions it 
corresponds to averaging over all ion trajectories. 

The fidelity of a two-qubit phase gate corresponding 
to the sequence S'{RG)'^ (i.e. phase gate with tt pulses) 
can be written as 



mm 



E 

a,/3=0 



(18) 



where a' = 1 — a and /?' = 1 — 0. Whenever SOap ^ 1 
and SQa/3 + ^^a'f)' ^^af) (achieved by applying the 
TT pulses), we obtain a significant improvement in the fi- 
delity of the phase gate using the n pulses, that is 
1 - J^' < 1 - JP. 
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B. Implementation on ions in microtraps 

Consider two ions confined in two separate harmonic 
potentials (microtraps) . with a logical-statc-sclcctivc and 
time-dependent external force acting on both ions. In 
the semiclassical approach (trajectories of the ions are 
considered to be classical) the Hamiltonian of the system 
is 



H{t) = J2 Ha0{t) \a)i{a\ 



(19) 



where 



Ha/sit) 



■ 



2m 2m 



1 



+ {s-Xa- d/2)F„{t) + (s' - 4 + dl2)F'i3{t) 
I 



+ T 



(20) 



where the two bare (empty) microtraps with a trapping 
frequency u) are separated by a distance rfo, m is the 
ion mass, Xa{t) and a;^(t) are coordinates (trajectories) 
of ions 1 and 2 corresponding to their internal states 
\a = 0, l)i and |/3 = 0, 1)2, Pa{i) a-nd p'p{t) are momenta 
of the ions, and d is the equilibrium distance between 
the ions (d > do) minimizing the total confining poten- 
tial (potential of microtraps -I- Coulomb repulsion). We 
denote i = q'^/Aneo, where q is the ion charge and en 
is the permittivity of vacmun. Finally, Fa{t) = aF{t) 
is a logical-state-selective (a = 0, 1) and time-dependent 
force which displaces the ion only when it is in its log- 
ical (internal) state \a = 1). The parameters s and s' 
are associated with the potential (introduced by the dis- 
placing forces F and F') which the ions experience when 
they are in their equilibrium positions {x = ±d/2) in 
the microtraps. 

The dynamics of the ion system during the G pulse are 
governed by the evolution operator 



U = Dexp 



Hit') dt' 



(21) 



where D denoted the Dyson time-ordering operator, and 
the Hamiltonian H{t) is defined in Eq. (19). The inte- 
gration is carried out over a time interval {to,t), where 
we will assume that to < 0, t > 0, |to|, \t\ ^ t, where t 
describes the diiration of the time interval on which the 
pushing force is applied. This assumption will become 
clear when we choose the time profile of the force later 
on in this section. Then, the corresponding evolution is 

U\a(3) = e'^-i'\a(3) , 



where 



Qa0 = -\f^H^p{t')dt'. 



(22) 



(23) 



The overall phase in Eq. (10) can be written (after a coor- 



dinate transformation ; 



Xa~d/2 and x'g x'^+d/2) 



as 



where 



-!? = 



ya/3 



'11 



+ < 



>00 ■ 



^Jto \d + x'i3-Xa 



■dt'. 



(24) 



(25) 



This means that the overall phase of the phase gate is 
determined only by the Coulomb interaction between the 
ions. 

Let us define an important parameter 



(26) 



2^3 ' 



which gives the ratio between the Coulomb repulsion en- 
ergy and the trapping potential. It is small (e 4C 1) when 
the traps (i.e. the ions as well) are far apart and it can be 
large (e ^ 1) when we move the traps close to each other. 
Here we consider the performance of the phase gate in the 
regime e <C 1. In Sec. V we extend the treatment to all 
values of e. 

In the regime e <C 1 an ion in the internal state |1) is 
displaced from its equilibrium position by 



x{t) 



F{t) 



(27) 



which also corresponds to the displacement of a single 

ion in a harmonic trapping potential (Appendix D in 
Ref. [6]). An ion in the internal state |0) is not dis- 
placed because the internal-state-selective force does not 
act on it. Thus, in the regime e <C 1 we can write for the 
displacement of an ion Xa{t) = ax{t), where a e {0, 1}. 
Then, it can be shown [2, 5] that the overall phase of the 
two-qubit phase gate in Eq. (10) is ■& f» 9 with 



(28) 



where we assumed a Gaussian time profile of the force, 
F{t) = Foe-(*/^)', introduced ^ = aFo/hcj, and de- 
noted a = ^h/moj. Then it follows from Eq. (27) that 

S(t) =^ae-(*/^)' . This means that the displacement x is 
measured in units a, where ^ is a dimensionless parame- 
ter. 

The time r for which the force needs to be applied is 
to be calculated from the condition ^ = n (i.e. when the 
phase gate is equivalent to the controlled-NOT gate). In 
practice, the internal-state-selective force F will be pro- 
duced in a non-resonant laser beam (dipole force) with 
F DC //A, where I is the laser intensity, and A is de- 
tuning from a driven atomic transition |1) <-> |aux). For 
convenience we will treat the case where the logical state 
|0) is not coupled to the laser, i.e. ensuring the state 
selectivity of the force. The extension to the case where 
both logical states experience a non-zero but diff'erent 
force is straightforward. 
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C. Cancellation of single-qubit rotations 

The TT-pulse method (spin-echo) improves the fideHty 

of the phase gate by cancclhng some single-qubit rota- 
tions which may be imprecise. Here, we propose two 
alternative ways to achieve the same effect without using 
any tt pulses. 

The phases Qap hi Eq. (3) have a structure [6] 

Qa0 = {a- P)ei ^ + (a - /3)2^2 + 0{a/d) , (29) 

where 61,62 are constants of order one (6*1,^2 ~ 1), we 
have typically a/d ^ 10~^, and we already know from 
the last paragraph in the previous section that ^ oc 
I/A. 

The second term in Eq. (29) is of order one because 
this term produces the two-particle phase (28), and the 
gate time is chosen to ensure this phase is tt or 7r/2. With 
^ ~ 1, which it will be seen later is a reasonable value, 
the first term is of order 10^. 

When the laser intensity fluctuates it causes the value 
of ^ to fluctuate as well. Let us assume, for example, 
1% laser intensity fluctuations. Then the first term in 
Eq. (29) produces an uncertainty of order ten, which is 
enough to completely spoil the gate. If this can be sup- 
pressed then the gate will work since the second term has 
a much smaller uncertainty. We note that the combina- 
tion 

601 + 610 = 2^2 + Oia^d^) (30) 

achieves the cancellation while preserving the $2 term 
which is needed for the gate. To obtain this cancella- 
tion, we can use two pulses with opposite signs for the 
61 term but the same sign for the 62 term. The 7r-pulse 
method discussed above achieves this by swapping the 
sign oi a — p. An alternative approach is to swap the 
sign of ^. Since ^ oc 1/A this can be done by using two 
successive pulses with opposite detuning A. This is our 
first suggestion. We thus replace the spin-echo pulse se- 
quence RGRG with the sequence GbGr, where Gb (Gr) 
is the G pulse with blue detuning (red detuning). 

We showed in Ref. [2] that the main contributions to 
9i are from the Coulomb energy and the light shift (AC 
Stark shift) produced by the laser field providing the 
dipole force. We there derived the condition for these 
to be equal and opposite (so called "sweet spot"). 

We next consider how to cancel them each separately. 
It is possible to make the contribution from the light shift 
zero by a judicious choice of laser polarization, such that 
the two internal states |0) and |1) experience equal light 
shifts but different forces [10, 11]. Only the Coulomb 
contribution then remains. This can be cancelled by 
applying the force in two successive pushes in opposite 
directions. The linear di term in Eq. (29) then cancels, 
while the 02 term remains. 

In all these methods, the cancellation takes place as 
long as the laser intensity is the same for the two closely- 
spaced pulses. The gate remains sensitive to intensity 



changes between one pulse and its partner in a given pair, 
but is much less sensitive to changes in laser intensity 
between one pulse-pair and another. 

Avoiding the spin-echo method is advantageous since 
the single-qubit 7r-pulses it requires may be slow or im- 
precise. However there are practical problems in all the 
methods. Switching of the detuning may be technically 
difficult since A can be of order GHz. Switching the force 
direction is not convenient in a travelling wave configura- 
tion where the force on each ion arises from the transverse 
profile of a single focused laser beam, however it can be 
done conveniently in a standing wave configuration by 
changing the relative phase of the two travelling waves 
forming the standing wave [10]. 

III. DECOHERENCE-FREE SUBSPACE 

Several advantages may be obtained by encoding log- 
ical qubits not in single ions but in pairs of ions, in the 
decoherence-free subspace (DFS) spanned by 

|0)^ = |01), |1)^^|10). (31) 

Here the subscript L refers to logical states in the sub- 
space, the kets without subscript refer to physical states 
of the trapped ions. This concept has been discussed in 
the context of trapped ions in Ref. [7]. The advantages 
arise from the fact that the two logical states evolve in 
the same way under any infiuence which causes the same 
phase change of |1) relative to |0) in both ions, therefore 
the logical information is completely protected from such 
infiuences. An example is a change in the magnetic field 
on both ions together, or a phase drift of the lasers used 
to implement the operations. 

G pulse G pulse 

ion qubIt logical qubit Ion qubit logical qubit 

(a) (b) 

FIG. 1: Two different ways of encoding logical qubits in 
pairs of ions qubits in the decoherence-free subspace (DFS). 
(a) Two pairs of ion qubits are arranged along a line to form 
a pair of logical qubits. We apply the G pulse only between 
the two central ions, (b) A more symmetric way of encoding 
logical qubits in the DFS. We apply the G pulse between the 
two pairs of ion qubits. 

We will next present two methods to implement the 
controUed-phase gate '~^Z in the DFS. The required evo- 
lution is 

|01)|01) ^ |01)|01) , 

|01)|10) ^ |01)|10) , 

|10)|01) ^ |10)|01) , ^ > 

mm ^ -iio)iio) . 



5 



First, consider two pairs of ions arranged along a line 
as in FIG. la. They may be all in separate microtraps, or 
with one or more in the same trap. All we need to do to 
achieve the desired gate is apply the pushing method to 
the central two ions. The presence of the other ions has a 
small influence on the distance through which the pushed 
ions are displaced, and their Coulomb interactions con- 
tribute to the total phase. However, this simply changes 
by a factor of order one the time required to achieve the 
desired phase tt, but it does not introduce sensitivity to 
thermal motion. We thus obtain the evolution 



|01) |01) 
|01)|10) 
|10)|01) 
|10)|10) 



|01) |01) , 
|01)|10) , 
|10) |01) , 
|10)|10). 



(33) 



To obtain the phase gate in Eq. (32), the Pauli oz oper- 
ation is then applied to the second logical qubit. 

It may be objected that the operation just described 
lacks symmetry in the way the controlling apparatus in- 
teracts with the ions in a given pair, and therefore it is 
against the spirit of the DFS concept. However, the joint 
state of the ions remains at all times in the DFS, and the 
mechanism of the pushing force (such as a far off-resonant 
dipole force) does not rely on the need for the very precise 
frequency matching which the DFS is designed to avoid. 
Therefore the need for very precise frequency matching 
(or magnetic field control) is still avoided. 

An alternative and more symmetric way to implement 
the gate in the DFS is shown in FIG. lb [13]. By pushing 
the ions as shown, before any single-qubit rotations are 
added the evolution is 



|01) |01) 
101) |10) 
|10)|01) 
|10)|10) 



ei(eoo+eii) |oi) |oi) , 

e'(^"i+«") |01) |10) , 
g.(eio+eoi) 1 10) |01) , 

e^^^^ii+Qof) |10) |10) , 



(34) 



where Qa{) are the same phases as in Eq. (3). Recalling 
Eq. (10), the overall phase for the evolution of the logical 
qubits is 

^ = (Gn + Goo) - 2(910 + eoi) + (600 + en) 

= 2(611 -eio-eoi + eoo). (35) 

Therefore the gate is obtained in a time half that required 
for single ions in the same conditions. 

The method of FIG. lb may be more difficult to re- 
alise in practice, owing to the two-dimensional rather 
than one-dimensional configuration of the ions. How- 
ever, it has the advantage that the cancellation of un- 
wanted phases in Eq. (29) takes place automatically with- 
out any need for tt pulses or switched detunings. This 
also means the "sweet spot" method discussed in Ref. [2] 
is not needed. 



IV. NON-NEIGHBOUR IONS AND TOFFOLI 
GATE 

A. Principle 

We now consider a system of three qubits. First we 

present an analysis of the general requirements along sim- 
ilar lines to that used for two qubits in Sec. II. Suppose 
a system of three interacting qubits evolves as 



|a/?7) 



la/37) 



(36) 



where |a/37) = |a;)i ® ® |c)3 denotes a joint logical 
state of three qubits with a, /3, 7 € {0,1}. Our aim is 
now to realize a three-qubit phase gate 

la/37) — ^ la/37), Q!/37 = 0, (37a) 
la/37) — ^ e^^\a^i), al3'y = l, (37b) 

that is the gate which produces an overall phase 'd if and 
only if all three qubits are in the logical state |1). The 
three-qubit phase gate can be obtained when we apply 
single-qubit rotations after the operation G, that is 



S = Si ® S2 ® S3 , 




(38) 


where 






5i = |0)i(0|e^^« + |l)i(l| 




(39a) 


= 10)2(01 e^^° + 11)2(11 




(39b) 


^3 = 10)3(01 e^^° + |l)3(l| 




(39c) 


The sequence SG will correspond to the three-qubit gate 


if the following equalities are satisfied 






^0 + -Bo + C'o + 6000 = 


0, 


(40a) 


Ai+Bo + Co + Qwo = 


0, 


(40b) 


Ao + Bi + Co + Ooio = 


0, 


(40c) 


Ao + Bo + Ci + Qooi = 


0, 


(40d) 


Ai + Si + Co + Olio = 


0, 


(40e) 


Ai + Bo + Ci + Old = 


0, 


(40f) 


Ao + Bi + Ci + 6011 = 


0, 


(40g) 



where the overall phase reads 

^ = Ai + Bi + Ci + 6111 . 



(41) 



In Eqs. (40) we have seven equations in six independent 

parameters (Ao,...,Ci), but the equations are not all 
linearly independent. By combining them we obtain 

6000 + 6110 = 6100 + 6010, (42a) 
6000 + 6101 = 6100 + 6001, (42b) 
6000 + 6011 = 6010 + 6001, (42c) 

These are conditions on the evolution (36) which will 
have to be satisfied if the Toffoli gate (the three-qubit 
phase gate, to be precise) is to be realised. 
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When we apply Eqs. (42), then Eqs. (40e)-(40g) can 
be expressed as hnear combinations of Eqs. (40a)-(40d). 
This means that the original system (40) is now reduced 
to four independent equations and six independent pa- 
rameters. Therefore, we can choose two parameters and 
calculate the rest. Let us pick Aq = Bq = — 6ooo/3. We 
then get Co = — 6ooo/3 and 



with 

Hal3j{t) 



Ai = -9100 + 26000/3, 
B\ = —0010 + 28000/3 , 
Ci = — 9ooi + 2©ooo/3 , 



(43a) 
(43b) 
(43c) 



+ 



{Xciij Xaj ) 
I 

\2d + xi;^-xa\ 



\d + x'^-x'f^\ ' 



(49) 



where the overall phase is 

= 6111 - 6100 - 6010 - 6001 + 26000 ■ (44) 

We can also introduce 7r-pulscs in order to cancel un- 
wanted terms. We replace the sequence SG by S'{RG)^, 
where .R = i?i (g) i?2 -R3 denotes a tt pulse applied on all 
three qubits. An analysis along similar lines to the above 
implies that the threc-qubit phase is now given by 

^' = 3(6000 + 6111) 

-(6100 + 6011 + 6010 + 6101 + 6001 + 6110) . 

(45) 

The single-qubit rotations S' = S'-^®S'<2^®S'o^ take the same 
form as in Eqs. (39), except that now we have = = 
Co = -(©000 + 6iii)/3, and 

A'l = -(6100 + 6011) + 2(6000 + 6iii)/3 (46a) 

B'^ = -(6010 + 6101) + 2(6000 + 6iii)/3 (46b) 

C'^ = -(6001 + 6110) + 2(6000 + 6iii)/3. (46c) 

There is again a set of conditions to be satisfied if the 
thrce-qubit phase gate with tt pulses is to be realized. 

They are 



where we introduced a notation ai = a, a2 = /3, as = 7 
u-«2 = a;^, Xas = a;", and we applied a 



:r" x'^ + d assuming that the ions are at sites x = —d, 
x = and x — d with e ^ 1. Then the three-particle 
dynamical phases Qa^-y read 



Qa0f 



Hal3-({t') dt' . 



(50) 



to 



The structure of the Hamiltonian in Eq. (49) allows us 
to write for the overall phase 



= (pill - 0100 - 9^010 - (pool + "^(pooo , 



(51) 



where 




d + x'o 



+ 



\2d + 



+■ 



d + 



dt'. 



(52) 



6110 + 6001 + 6000 + 6111 

6101 + 6010 + 6000 + Qui 



6 



100 



6 



oil 



6, 



010 



6011 + 6100 + 6000 + 6111 



6101 , 
(47a) 

6100 + Qoil + ©001 + ©110 , 

(47b) 

6010 + ©101 + 6001 + 6110 
(47c) 

Thus, the general concept of the three-qubit phase gate 
(which is equivalent to the Toffoli gate) seems to work in 
both cases. The question is now whether we are able to 
satisfy Eqs. (42) or (47), respectively, on the system of 
ions in microtraps. We discuss this issue in the following 
section. 

B. Implementation on ions in microtraps 

The Hamiltonian for the system of three ions in three 
separate microtraps is 
1 

H{t)= J2 i?a/37WI")l(«l®l/?)2(/3|® 17)3(71, (48) 
a,/3,7=0 



It follows from Eq. (51) that the overall phase is deter- 
mined again only by the interaction phases (pafj-,- rather 
than by the total phases Qafij- However, in contrast 
to Eq. (25) this time we have to consider three pairs of 
mutual Coulomb interactions. 

We need to verify whether Eqs. (42) can be satisfied 
for the system of ions in microtraps. In the adiabatic 
approximation {cut ^ 1) and for e -C 1 we can write the 
trajectories of the ions in the form 



x'0 



Xa + S , 
x'0+6', 

x'' + 5", 



(53a) 
(53b) 
(53c) 



where x^ denotes the state-selective displacement of the 
ion (Eq. (27)), and 5 corresponds to thermal oscillations 
of the ion in the microtrap. The adiabatic regime means 
that the state-selective pushing force F is applied on a 
time scale much longer than an oscillation period of the 
ions in microtraps. We assume first of all that state- 
dependent forces of the same magnitude and the same 
pulse time r are applied to all three ions. Then, the 
analysis of Eq. (52), using the methods of [2, 6], shows 
that the three-particle phases take the form 



(a-/3)2 + l(a-7)2 + (/J-7)2 



(54) 
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where 6 is defined in Eq. (28), ct,i3,'y G {0,1}, and we 
dropped terms of order 0{a^/(f) because we have typi- 
cally a/d ~ 10-2 - 10-3. 

When we apply the gate pulse G on three ions for a 
time T, followed by the single-qubit rotations S, we ob- 
tain 

^ e*^|110), 
-e^^/«|101), 
^ e'^lOll), 
^ e^i7^/8|lll), 

which corresponds to '^^12(6') '^^23(6') '^^'13(^/8), where 
'"PijiO) is the phase gate in Eq. (1) between ions i and j. 
This means that the single pulse G on all three ions simul- 
taneously (followed by single-qubit rotations S) produces 
three two-qubit controlled-phase gates. The two-qubit 
phase of the gate between the outermost pair of ions is 
6'/8, rather than 6, because the separation of those ions 
is 2d and the two-qubit gate phase is proportional to the 
inverse cube of the separation. 

To obtain the simple thrce-qubit phase gate '-^'^Z from 
Eq. (55) we apply a network as in FIG. 2. A sequence 
containing two controlled-NOT and two controlled-phase 
gates on adjacent ions causes the overall evolution to be 



(56) 



which is a three-qubit phase gate Z for 9 = A-k (or an 
integer multiple of 47r). The controlled-NOT gates can be 
obtained from the two-qubit phase gate plus Hadamard 
rotations (or 7r/2 pulses). 

Note that when = Air, the operation (55) is 
'^Pi3(7r/2). In other words, we have simply a two-qubit 
phase gate between the outermost pair of ions. In this 
case the force on the central ion is redundant, i.e. the 
same result would be obtained with laser beams pushing 
only the two outermost ions. This also shows that the 
network of FIG. 2 is a standard decomposition of ^^Z as 
presented in Ref. [12]. 

The implementation of the phase gate including a spin- 
echo by TT pulses is also completed by the same network 
(FIG. 3). 

We note that in principle, two-qubit gates between 
pairs of ions even further apart could also be implemented 
simply be pushing the relevant pair of ions, without the 
need first to make them adjacent. 



C. Speed of the gate 

Next we estimate the speed of three-qubit phase gates 
with ions in microscopic traps. We will consider single- 
qubit gates (tt pulses, rotations) to be fast compared to 



|000) - 
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FIG. 2: Three-qubit phase gate. A SG pulse sequence fol- 
lowed by a quantum-gate network produces a three-qubit 
phase gate when 9 = An. The network corresponds to the 
transformation in Eq. (56). 
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FIG. 3: Three-qubit phase gate with n pulses. A S'{RG)^ 
pulse sequence, followed the same network as in FIG. 2, pro- 
duces a three-qubit phase gate with tt pulses, i.e. a gate with 
a better performance in terms of fidelity. Note that each G 
gate in this network is obtained by applying the pushing force 
for half as long as for G in FIG. 2. 



multi-qubit gates. This is typically the case for trapped 
ion qubits. Therefore we only need to add together the 
times for two-qubit and three-qubit logic gates on the 
ions. We assume that these are all produced by the state- 
selective pushing method under discussion. 

There are two different assumptions about the speed of 
the operations which are both physically sensible. First, 
we may assume that all operations involve pushing forces 
of the same magnitude. In this case the time required to 
produce the two-qubit phase gate between the outermost 
ions is eight times longer than the time to produce the 
same gate between adjacent ions. On the other hand, 
we might assume that the forces are always adjusted in 
magnitude so that the gate duration is limited by the 
adiabatic criterion ujt ^ 1. In this case the gate between 
outermost ions has the same duration as the same gate 
between adjacent ions, but requires larger forces. 

With the former assumption (forces of given magni- 
tude), the total time required for the whole network, in- 
cluding G and the two-qubit gates, is 



Ttotal ^4t + T+ — 



3r 
2" 



(57) 



where r is the time required for a '-'Z gate (i.e. 9 = tt) 
between adjacent ions. The duration of the G pulse is 4t 
because we need 6 = An. The duration of the final gate 
is t/2 since exp{i<p) = exp[i{(t> + 2n}]. 

With the other assumption (all gates of duration t), 
the total time is 5t. 

The same conclusions apply to the three-qubit phase 
gate with tt pulses in FIG. 3, except that now each of the 
two G pulses is applied for a time 2r (instead of At), and 
the fidelity of the gate is higher. 
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V. TWO-QUBIT PHASE GATE IN GENERAL 

In Ref. [2] we discussed the dynamics of a two-qubit 

phase gate with ions in microtraps in the regime e ^ 
1, that is when the two microtraps were sufficiently far 
apart. In what follows we will re-examine the dynamics 
of the ion system outside this regime, i.e. e ^ 1. We 
will use a more general approach compared to the one in 
Ref. [2]. The aim is to understand the phase gate for all 
values of e, and especially for the case of two ions in the 
same trap, where e = 2. 

To this end it is useful to write the Hamiltonian in 
Eq. (20) in a new coordinate system 



Ra0 = {xa + a;^)/2 , 



(58a) 
(58b) 



where R^/s is a ccntrc-of-mass coordinate of the ion sys- 
tem, raj3 is a relative excursion of the ions from their 
equilibrium positions, and a,/3 G {0,1} correspond to 
logical states of ion qubits |0), Then, the Hamilto- 
nian splits into a centre-of-mass and a relative-motion 
part 



Ha/3 = H, 



Q/3 

2M 



a/3 



2n 



Va0 ) 



(59) 



where we have introduced Pa^ = MRa/s, Pap = /tw^a/s, 
M = 2to, and ji = to/2. The potential-energy terms of 
the two ions in a joint internal state \aP) are 



1 



= -Mu,'Rl0-Ra0{Fa+F'^) + {sFa + s'F'^), 



(60a) 



1 



= ^l^u:\rap + Mf - -rap{F'^ - F^) + , 

(60b) 

where we denoted Ac? = d — do (i.e. the difference be- 
tween the separation of ions loaded in the traps and the 
separation of the centres of bare traps), and we assumed 
\d + ra/3\ = d + raf3 since d » {rapl > 0. 

The ion separation d can be calculated from the equi- 
librium when no force is applied, that is 



(i) When 77 ^ (which corresponds to — > 00, i.e. mi- 
crotraps are far apart), we get e = q'^ /(Treomuj^dQ). This 
means that e <C 1, and d = dg. (ii) On the other hand, 
when 7] ^ 00 (corresponding to do 0, i.e. the two sep- 
arate microtraps overlap completely and we end up with 
two ions in a single linear trap), we get e = 2, which is 
also the maximum value of e. 

Substituting e = 2 into Eq. (26) gives 



d = 



1/3 



(63) 



which is the standard expression of the separation of two 
ions in the same harmonic trap. We refer to this case as 
the case of a "linear trap" , since typically one chooses a 
linear geometry in the Paul trap in order to avoid micro- 
motion when more than one ion is trapped. Note that 
whereas in separate microtraps we could adjust the trap- 
ping frequency and ion separation independently, in a 
linear trap (i.e. two ions in the same trap) they are mu- 
tually dependent. 



A. Anharmonicity 

In the discussion about ions in two separate microtraps 
in Sec. II, the results were valid in the regime e <C 1. Next 
we examine the meaning of this condition, and show how 
it can be substantially relaxed. This is necessary in order 
to allow the gate to be implemented for two ions in the 
same trap, since then e = 2. 

We can expand the Coulomb term in Eq. (60b) in a 
Taylor series owing to \ra0\ <C d, giving 



a/3 



-22 
2^'^ ^a/3 



1- 



r-a/3 



e + 1 d 



--r^p{F'p-F^), 



(64) 



where we denoted u) — + e, used the definition of e 
in Eq. (26), applied Eq. (61), and dropped a constant 
term Vq = fJ^uj^ (Ad)'^ / 2 + £/d. This means that in the 
regime e ~ 1 the effective trapping potential becomes 
anharmonic. We need to find out how the anharmonicity 
affects the centre-of-mass and the relative motion of the 
ions. 



dV 
dr 



IiuPM - — = , 
d^ 



givmg 



Ad: 



sinh^ \ -\n 
3 I 6 



'n + l + ^/ri(riT2) 



(61) 



(62) 



where 77 = £/[muP{do/'i)^\. Using Eq. (62) and substi- 
tuting d = do + Ad into the expression of e in Eq. (26), 
we find that e is a function of a single variable 77. The 
function is monotonic and there are two limiting cases. 



B. Dynamics 

Eqs. (60) and (64) lead to the equations of motion 

Rc0+u)'^Ra0 = Fan/M, (65a) 
+ ^1 - 7-^^^^ = , (65b) 

where F^f^ = F^ + F'^, = (F;, - F^)/2 with a,0 e 
{0, 1}. In the following we assume for convenience that 



9 



the statc-sclcctive-force Fa acts only when the ion is in 
its internal state \a ~ 1), so that F^p = {a + P)F, f^p = 
{13 - a)F/2, where F = Fi = F[. Using Appendix D in 
Ref. [6] the solution of Eq. (65a) in the adiabatic regime 
(wr » 1) is 



i?a/3W«i?a/3(i) + A(t), 



with 



Ma;2 



x{t) 



2Er 
Mw2 



cos(a;i + ^tRj , 



(66) 

(67a) 
(67b) 



where we used Eq. (27), and Eji {iPr) is the oscillation 
energy (initial motional phase) of the centre-of-mass mo- 
tion of the ions. 

There is no analytical solution of Eq. (65b). However, 
owing to the fact that Ivapl "C d the anharmonicity is 
small even for e 1. Therefore, we can estimate the 
solution in the adiabatic approximation {ojt ^ 1) as 



(68) 



with 



(/3-a) 



3e fa^jt) 
2(e+l) nuj'^d 



e + 1 



l2E, 



Scxffit) = \ —::^COs{u)t + tpr) + Cafs{t): 



(69a) 
(69b) 



where fap arises from the condition dVaf}/ drap = 0, Ej. 
and refer to the oscillation energy and an initial mo- 
tional phase of the relative motion of the ions, and Ca/3 
is to be discovered. 

When only one ion is in the state |1), that ion is 
displaced by the quantity x{t) imposed by the force in 
Eq. (27). It follows from Eq. (69a) that when e < 1 
this is also the change in the relative separation, i.e. the 
other ion hardly moves besides random thermal oscilla- 
tions. However, in the regime e ~ 1 the ions are close 
enough to be strongly coupled. Then the other ion moves 
as well and the net result is that the relative displacement 
is reduced by a factor 1 -|- e. 

The main consequence of anharmonicity is the extra 
contribution C,afj{t) to the relative motion of the ions in 
Eq. (69b). We no longer have the special property of a 
harmonic potential that the average excursion under free 
oscillatory motion is zero (when averaged over a large or 
integral number of oscillations). Instead the average is 
(Ca/s)- We do not have an analytical solution and so we 
must estimate Ca/3(*)- As an estimate we use the average 
of the two turning points of the motion (at given energy) 
in the anharmonic well, i.e. = (^as + ''afl)/^' 'where 



' a/3 



and r"^_j are the two roots (closest to the equilibrium 



point Tq,^) of the equation 



V^p{r^0)=Er + V^p{f^0). 



(70) 



where V^^ is the anharmonic potential defined in 
Eq. (64), and Er is thermal energy of the relative mo- 
tion of the ions in Eq. (69b). We expect this estimate 

to be good up to a numerical factor of order one. By 
obtaining the two roots and expanding C,ai3 as a Taylor 
series in Er, we find 



WW 



Er 



e+1 jjLUj'^d 



1 + 



6e Japit) 



e+1 /UfD^d 



(71) 



Note that the influence of Ca/3 on the phases appearing 
in the quantum logic gate arises only from that part of 
which depends on a and /?, the rest will merely add 
a global phase. A numerical solution of the equations of 
motion confirmed the main points of this estimate. It 
was found that (Ca/s) typically underestimates {Sap) by 
a factor approximately six. 

In order to calculate the dynamical phases in Eq. (23) 
it is crucial to include both the centre-of-mass motion and 
the relative motion, in particular because the two-qubit 
phase is found to arise from a difference between the two 
contributions, associated with the fact that uj ^ uj. By 
substituting Eqs. (66) and (68) into the Hamiltonian (59) 
we determine the dynamic phases in Eq. (23) and hence 
obtain the overall phase in Eq. (10). The result is 



e + 1 



1 



1 2 + e 
{ujtY 1 + e 



+ 0(aV(i') 



(72) 



where 6 is given by Eq. (28). 

In the limit when e <C 1 and lot ^ 1 we obtain Ki 9 
as before. The main effect of bringing the traps together 
is to introduce the factor l/(e + 1). The (small) term in 
l/(a;T)2 in Eq. (72) is a surprise since it survives even 
when e — > 0, and therefore it should have been present 
in the earlier calculations. It arises from the fact that 
the kinetic energy of the vibrational motion depends on 
where the ion is in the harmonic well, and we now allow 
for the fact that Co ^ to in all kinetic energy terms whereas 
previously we did not. 



C. Phase gate in a linear trap 

Lot {)l be the overall two-qubit phase from Eq. (72) in 
the case of two ions in the same linear ion trap (e = 2). 
Then 



(73) 



If there is a sufficient force available (represented by ^i), 
then the gate time is limited by the criteria that (i) the 
ions are not displaced too far [x/d « aS^h/d <C 1), and 
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that (ii) the adiabatic conditions foUowing from Eqs. (66) 
and (68), that is Cjtl ^ 1 and ujtl ^ 1, hold. TypicaUy, 
the latter is the limiting condition, and it is sufficient 
to choose LOT = ojtl > 5 because to be precise what we 



require in practice is e 



rather than lut 3> 1. It 



follows that for a given value of ujt to achieve a desired 
overall phase the pushing force is smaller when the ions 
are in a linear trap (e = 2) compared to when they are 
in two separate microtraps (e < 2). The ratio is 



(74) 



This can be useful to reduce photon scattering (see be- 
low). 

The infidelity of a two-qubit phase gate for ions in two 
separate microtraps with e <i; 1 was discussed in detail in 
Ref. [2]. The dynamics when we consider general e was 
presented in Sec. V. This permits us to make a direct 
comparison between the case when the ions are in sepa- 
rate traps (e <C 1) and in the same trap (e = 2), and thus 
carry over the results of the previous study in Ref. [2] to 
the present one. 

We assume the state-selective pushing force is the op- 
tical dipole force produced by an intense laser beam (see 
the last paragraph in Sec. II B) in either the travelling- 
wave (TW) or the standing-wave (SW) configuration. 
There are two main contributions to the total infidelity 
of the gate: (i) thermal averaging and (ii) photon scat- 
tering. 

Two types of thermal averaging take place. These are 

the averaging over the dynamic phases Oq,^ (2?), and 
the averaging over the spatial dependence of the pushing 
force (S). We show in Appendix A that the contribu- 
tion 2? is of a similar order of magnitude for ions in sep- 
arate traps and in the same trap. We showed in Ref. [2] 
that the contribution S dominates "D for ions in separate 
traps when the waist of the laser beam is w < d (TW), 
or the laser wavelength is A < (SW). It follows that the 
same conditions cause <S to dominate V for ions in the 
same trap. 

The averaging over the spatial profile of the force (<S) 

is associated with the thermal motion of either ion with 
respect to the laser field illuminating it. We assume for 
the sake of argument that for the thermal oscillations of 
two ions in the same trap, the centre-of-mass mode (uj) 
and the breathing mode (wVS) have the same tempera- 
ture. Part of the random phase associated with thermal 
motion may be cancelled by a common-mode rejection, 
but to be cautious here we will ignore that possibility 
and therefore we overestimate the infidelity. In this case 
the infidelity associated with process S is approximately 
the same as calculated in Ref. [2] for given values of trap 
frequency uj, temperature T and laser parameters. We 
assume the spin-echo sequence (tt pulses) is used, so we 
have two gates each with i? = 7r/2, and we add their two 
infidelities. For travelling wave excitation with Gaussian 
laser beams (each ion illuminated by a separate beam) 



the resulting infidelity is [8] 



TW 



3 \ huj 
2 /STrfcsT 



( 2a;o w 
wJ \ w 2xo 



(75) 



where Q{y) = Uy'^ - 64y'^ + 89- 34/^^ + 1/y'^, oj is the 
trapping frequency, a = ^/h/rmv, w is the size of the 
waist of the laser beam, and xq is the position of the 
ions in the profile of the laser beam. Placing the ions at 
Xq =w/2 offers a useful improvement in the performance 
of the gate. 

For standing wave excitation (i.e. the pair of ions po- 
sitioned in an optical standing wave along z) we have [9] 



128 



1 — exp 



(76) 



where ka = (47r/A) sin(Q:/2), a is the the angle between 
two laser beams forming the standing- wave field, zq is the 
ion position in this field, and we assumed that the ions 
are placed at kaZo = 7r/4 which minimizes the infidelity 
and simplifies its expression. 

The expression for the number of scattered photons 
N during the gate operation in a linear trap is de- 
rived in the same way as the result for ions in sep- 
arate traps in Ref. [2], except that now we consider 
tl = 3^/77/(0)^1 \/2), where we assume the spin-echo se- 
quence (hence t) — tt/2) and we add the numbers of scat- 
tered photons in the pair of gates. Treating the electronic 
transition giving rise to the dipole force by a two-level 
atom model, the result for travelling-wave excitation is 



TW 



where P is the laser power, and 



C" 



7r*c m 
2^2 A3 



(77) 



(78) 



where c is the speed of light, m is the ion mass, and A is 
the laser wavelength. 

For standing-wave excitation we obtain 



io' 1 

P COs'^{kaZo) ' 



where [14] 



C" 



TT c m 1 



(79) 



(80) 



4^2 A sin2(a/2) " 
The total infidelity of the phase gate can be expressed 



as 



tot 



(81) 
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assuming V,N <^ 1. 

FIG. 4 and FIG. 5 give the total infidelity for some 
example parameter choices, with atomic properties ap- 
propriate to the ^°Ca+ ion. 
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FIG. 4: The total infidelity of a two-qubit phase gate for two 
ions in a linear trap as a function of the trap frequency for 
various parameter values for travelling-wave excitation. The 
solid lines are for w = 4 /im and P = 10 mW, the dashed lines 
are for w — 2/im and P — 100 mW. The thick lines are for 
the temperature equal to the Doppler temperature (Tdopp = 
538 /xK) associated with the dipole-allowed 451/2 ^ 4Pi/2 
transition at A = 397 nm in *''Ca+. The thin lines are for 
T = huj/{kB In 2). The gate time is tl = S/o;. 
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FIG. 5: The same as in FIG. 4 but for standing-wave excita- 
tion, with a. — -R /2 and k^ZQ = 7r/4. 



VI. CONCLUSION 

In this paper we offered further quantum-gate methods 
using the state-selective displacement of trapped atomic 



ions [1], and in this sense this paper complements our 
previous publication [2]. 

We proposed methods to suppress unwanted effects in 
the evolution, as an alternative to the spin-echo method. 
We also discussed the application of the gate to logical 
qubits encoded in pairs of ions (rather than in single ions) 
in a decoherence-free subspace. 

We showed that it is possible to realize two-qubit 
gates between non-neighbouring ions without the need 
to swap the ions around or move quantum information 
between them. This produced an implementation of the 
Toffoli (three-qubit controUcd-NOT) gate with three ions 
in three separate microtraps, where the three-qubit gate 
is five to eight times slower than its two-qubit counter- 
part. 

We also showed that the original state-selective dis- 
placement method for a two-qubit phase gate with ions 
in separate microtraps can be extended to the case with 
two ions in the same trap. We analyzed the anharmonic- 
ity of the effective trapping potential, arising from the 
Coulomb interaction potential between ions, and found 
that its effects were small even in the limit of two ions in 
the same trap. Therefore the gate retains its attractive 
features (a good combination of speed and robustness), 
whatever the separation of the ion traps. 
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APPENDIX A 

We showed in Ref. [6] that in the regime e <C 1 a con- 
tribution from thermal averaging to the total infidelity 
of a two-qubit phase gate (defined in Eq. (81)) is given 
by Eqs. (75) and (76). In this contribution thermal av- 
eraging over the size of the state-selective force (i.e. over 
the spatial profile of the force) dominates averaging over 
the dynamic phase 9q/3 (i.e. over a discrepancy between 
random values of Qai3 in a G pulse and their determin- 
istic values in single-qubit rotations S), where the latter 
reads [15] 

-.«-(^)^(^)'- <-) 

In the regime of interest in this paper (e = 2) we need 
to reconsider this conclusion on the grounds of anhar- 
monic effects present in the system. In particular, we 
need to calculate the contribution to the infidelity from 
thermal averaging over the dynamic phases for e — 2 and 
compare the result to the contribution from averaging 
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over the force profile, which is the same for any value 
of e. 

Using a precise analytical expression of the overall 
phase 1? in the regime e = 2 we calculate that the in- 
fidelity from averaging over the dynamic phases is 




When wc use the phase condition = tt, and choose 
u>T = 5, we obtain 




Thus wc find that the contribution arising from an- 
harmonicity scales as (a/d)^, and therefore it docs not 
dominate the infidelity. Comparing Eqs. (Al) and (A3), 
the expressions have the same functional form and a sim- 
ilar size. Both are small compared to Eqs. (75) and (76) 
when the laser beam waist or the standing wave period 
are small compared to the ion separation d. 
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